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Abstract 

We describe a new method to determine the position of the atomic cloud during its interaction 
with the microwave field in the cavity of a fountain clock. The positional information is extracted 
from the spectrum of the \F = 3, mj? = 0) to \F = 4, mp = — 1) hyperfine transition, which shows a 
position dependent asymmetry when the magnetic C-field is tilted by a few degrees with respect to 
the cavity axis. Analysis of this spectral asymmetry provides the horizontal center-of-mass position 
for the ensemble of atoms contributing to frequency measurements. With an uncertainty on the 
order of 0.1 mm, the obtained information is useful for putting limits on the systematic uncertainty 
due to distributed cavity phase gradients. The validity of the new method is demonstrated through 
experimental evidence. 
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I. INTRODUCTION 



Caesium fountain clocks provide the best available implementation of the SI second with a 
relative uncertainty below 10~ 15 [1]. As such they are crucial to the accuracy of International 
Atomic Time (TAI). 

The basic experimental setup of a fountain clock is well established: A large number of 
atoms is laser-cooled to /iK temperatures. The atomic cloud is then launched vertically, 
passing through a cavity during its upwards motion. Here it is exposed to a microwave field 
that drives the clock transition and puts the atoms into a quantum mechanical superposition 
state. Gravity causes the atoms to slow down and fall back, where they cross the same cavity 
again and experience a second microwave pulse. This sequence implements Ramsey's method 
of separated oscillatory fields with a single cavity, avoiding any frequency shift resulting from 
a global phase difference between the two field regions. However, the atoms generally traverse 
the cavity at different horizontal positions during the upwards and the downwards passage, 
such that position dependent phase variations of the microwave field cause a frequency 
bias. This gives rise to one of the leading contributions to the uncertainty budget of several 



atomic fountain c 
phase gradients 



ocks An improved understanding of the actual microwave fields and 
has recently led to a renewed investigation of these effects 4|-[6| and 
it has been shown that under normal operating conditions the largest frequency shifts can 
be expected to occur due to phase gradients oriented horizontally across the cavity. Such 
gradients result from imbalances in the microwave feeds or from inhomogeneities of the 
electrical conductivity of the cavity material. 

An assessment of the resulting contribution to the frequency uncertainty requires infor- 
mation on the horizontal position of the cloud during its upwards and downwards cavity 
passage. So far, no method for an in-situ measurement has been available. Estimates from 
images of the atom cloud before the launch, or possibly during the passages through the de- 
tection zone, carry a particular difficulty: Information is needed on the distribution of those 
atoms that are detected at the end of the sequence, and this can differ significantly from 
the full distribution of all launched atoms. The full atomic trajectories and their interaction 
with the known obstructions inside the fountain can be taken into account by simulations, 
but the results strongly depend on the assumed cloud parameters. 

For the majority of cavity designs in use, it is possible to deliberately apply a phase 



gradient in the cavity by asymmetric feeding. This allows an adjustment of the launch 
direction to effectively cancel the shift of the cloud position between passages for one coor- 
dinate direction 4, 5|. Improved cavity designs would permit this for both directions j^], 
but have not been implemented in any fountain yet. Furthermore, this method requires 
manual adjustment of the fountain alignment and considerable measurement time, making 
it unsuitable for quick confirmation measurements or even on-line monitoring of the cloud 
position. The cavity design used in the PTB fountains {t| favors stability of the microwave 
feeding over adjustability and does not provide this option at all. 

In this paper, we present a new method to find the center-of-mass position during the 
upwards and the downwards cavity passages for the same fraction of atoms in the cloud 
that also contributes to the frequency measurement. The method is based on position 
dependent asymmetries in the transition probability spectra of the |-F = 3, = 0) to 
\F = 4, rriF = — 1) transition. Its use in 6] has helped reduce the uncertainty of the fre- 
quency bias from the cavity phase distribution in PTB's fountain CSF2 to 1.33 x 10~ 16 and 
the total systematic uncertainty to Ub — 4.1 x 10~ 16 . 

Section [Til of this paper aims to provide an illustrative explanation for the asymmetry of 
the measured spectra and its relation to the magnetic field in the cavity. In section [TTTI these 
considerations are developed into an approximated analytical model. This provides a simple 
way to determine the orientation of the magnetic field and to extract the center-of-mass 
position for the ensemble of atoms contributing to frequency measurements. The results of 
these calculations are consistent with experimental data, as presented in section IIVI The 
uncertainty of the position determination is discussed in section |V] 



II. ATOM-FIELD INTERACTION FOR Am = ±1 TRANSITIONS 

In the fountain CSF2, the atoms are state-selected for the \F = 3,mp = 0) state by a 
combination of microwave interaction and light pressure before entering the cavity [§]. A 
nominally vertical magnetic field, usually referred to as the "C-field" , makes the transitions 
to different Zeeman sublevels of the \F = 4) hyperfine state individually addressable. We 
will use the notation [a to b] to refer to transitions \F = 3, m F = a) to \F = 4, m F = b). 
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A. Field distribution in a fountain clock cavity 



Fountain clocks use cylindrical cavities designed such that the TE n mode is resonant at 
the frequency of the [0 to 0] clock transition. The atoms enter and leave the cavity through 
circular apertures in the top and bottom endcaps. The diameter of these apertures is chosen 
below the wavelength of the microwave signals such that the field distribution remains 
confined. The addition of cutoff tubes of the same diameter further reduces microwave 
leakage, which might otherwise lead to a bias of the measured clock frequency (9), [uj] . 

The actual fields in the cavity, including the effects of endcap holes and wall losses, 
can be calculated using finite element methods 

as 

. For an analytic description, we will 
approximate the field distribution by that in a perfectly conducting cavity without holes. For 
a cavity with height d centered around the origin, this is described in cylindrical coordinates 
by the equations 
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The Bessel functions Jo an d J\ describe the field dependence on the radial coordinate p, and 
7 = Xo/R, where xq is the first zero of J\ and R is the cavity radius. In the PTB fountains 



We will take the equations to describe 
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CSF1 and CSF2, R = 24.2 mm and d = 28.12 mm 
the instantaneous field distribution at t — 0, and they have been normalized such that at 
p = the integral of H z over the vertical coordinate z is 

'd/2 
-d/2 

The approximation reproduces the key features of the field distribution very well: When 
traveling vertically through the cavity, the sign of the horizontal field amplitudes H x and 
H y reverses at the center plane. Moving radially, H x and H y are zero on the cavity axis 
and increase with radial position p. In the inner region that is accessible to atoms traveling 
through the endcap apertures, they are generally small compared to the vertical amplitude 
H z , except at the top and bottom, where H z becomes small. Radially, H z has its largest 
value on the cavity axis and falls off slightly with increasing p. 

An atom is assumed to enter the cavity at —t c /2 and leave it again at +t c /2, with 
t c = 10.6 ms for CSF2. During the traversal, the motion of the atoms translates the spatial 



distribution of the field into a temporal pulse. Due to the extremely narrow linewidth of 
the clock transition, the amplitude and phase of this pulse determine the features of the 
transition probability spectrum. 



B. Angle — phase equivalence 

The most significant phase effect results from the local direction of the rf field. An 
oscillating horizontal field component oriented at an angle to the x-axis and with a peak 
amplitude Hq, can be written as a vector in xy-coordinates: 

mm = ?i ( -^A + ^ , . (4) 

2 I- sin[w rf t- /3] I 2 I sin[u; rf t + 0\ 
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Physically, this treats the linearly polarized field as the superposition of two opposite circular 
polarization contributions H cw and H ccw . Each of these carries the angular momentum along 
the nominal C-field axis that is required to drive one of the Am = ±1 transitions. We have 
chosen to investigate the [0 to -1] transition, which has the same initial \F = 3,m^ = 0) 
state used in clock operation and interacts with the clockwise circular contribution H cw 



(e.g., see 111]). In the relevant term, the angle (3 appears as a phase retardation. In essence, 
atoms experience any change in the orientation of the linearly polarized rf field as a change 
in its phase. 

The H ccw contribution, for which (3 advances the phase instead, can only excite the [0 to 1] 
transition. The Zeeman splitting is usually large enough that the off-resonant excitation of 
this transition does not add significantly to the transition probability while measuring the 
spectrum of the [0 to -1] transition. 



C. Shape of Am = ±1 spectra 

We now investigate the changes in the direction of the rf field experienced by an atom 
crossing the cavity. Let us assume that the reversal of the horizontal rf component between 
the top and the bottom of the cavity occurs in a smooth clockwise rotation. Note that this 
rotation is slow compared to that of H cw , and the field is best described as having a linear 
polarization with a slowly varying direction f3[t] = —irt/t c . This slow rotation causes a 
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gradual phase shift of the clockwise contribution driving the [0 to -1] transition, which can 
be expressed as a shift of the effective rf frequency by u rot = d(3/dt = iz/t c . The appearance 
of frequency shifts for the Am = ±1 transitions when the direction of the microwave field 
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changes relative to that of the static magnetic field is known as the Millman effect 

However, the horizontal field component normally does not rotate during the cavity 
traversal, but simply decreases to zero amplitude at the center plane and then reappears 
oriented in the opposite direction. This behavior can be modeled by taking the previously 
assumed linearly polarized field with its slow clockwise rotation and adding a second lin- 
early polarized field with a slow counter-clockwise rotation. This additional field then has 
its own fast-rotating, circular contributions, with only the clockwise contribution driving the 
[0 to -1] transition. The opposite direction of the slow rotation reverses the phase shift with 
time, such that the effective frequency of the second field has an opposite shift of — u rot . 

The presence of these two components in the rf field affecting the [0 to -1] transition 
suggests a single-passage transition probability spectrum with two peaks near ±u; ro t rela- 
tive to the unperturbed resonant frequency of the atoms. This is quite close to what is 
observed [ijj]. In CSF2, where u Iot = 7r/t c = 2n x 47.2 Hz, the two main lobes appear at 
±27r x 64.6 Hz due to interference effects that are not covered in this simple approach. 



D. Microwave pulse in the presence of a C-field tilt 

If the C-field is tilted by an angle a relative to the vertical z-axis, the rf components that 
are orthogonal to the field depend not only on H x and H y , but also on H z . We choose the 
coordinate system such that the C-field tilt is towards the negative x-direction, as shown 
in Fig. [TJ Additionally, we define a coordinate system (x',y,z') that is tilted around the 
common y-axis by the tilt angle a, such that the [0 to -1] transition is now driven by the 
unchanged H y component and H x i = H x cos a+H z sin a. Since H z is generally much stronger 
than H x , it can have a significant effect even for small tilt angles. 

This is illustrated in Fig. HJ In the cavity's center plane at z = 0, the horizontal com- 
ponents H x and H y are zero. The field is then given by the projected H z component and 
therefore oriented towards +x' everywhere. Near the top and bottom of the cavity the 
horizontal components dominate over the weak H z field and the radial symmetry is mostly 
preserved: The field points inwards at the bottom of the cavity and outwards at the top. 
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FIG. 1. Mixing of cavity rf field components: The tilt a of the C-field defines a new coordinates 
system (x',y,z') (note the shared y-axis). In this system, the field component H z r, consisting of 
the projected components flyw (a) and H z im (b), drives the [0 to 0] transition (a). [0 to -1] is 
driven by the unchanged H y component and H x /, consisting of the projections H x >[ z ] (a) and H x i^ 
(b). 

The transitions between the top, center and bottom regions are smooth, which generally 
causes the direction of the linearly polarized field experienced by an atom to rotate during 
a cavity passage. For an upwards passage at y < 0, the rotation is clockwise, resulting 
in a positive shift of the effective microwave frequency driving the [0 to -1] transition, as 
discussed previously. While the rotation in most cases does not occur at a constant angular 
velocity, this still breaks the symmetry of the spectrum around the resonant frequency and 
tends to lead to an enhancement of the lobe near — u TOt . For a traversal at y > the rotation 
occurs in a counter-clockwise direction instead, leading to an enhanced lobe near +w rot . 

The reversal occurs without rotation only for a passage at y — 0, where the projected H z 
component is parallel to the radial component, such that its only effect is to shift the point 
of zero field amplitude away from z = 0. 

III. ANALYTICAL MODEL 

In the following, we will develop a model based on a first-order quantum mechanical 
treatment and the analytical closed-cavity fields. This shows a simple way to quantify the 
asymmetry of the [0 to -1] spectra, which allows us to extract the center-of-mass position 
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FIG. 2. A tilt of the C-field results in a position dependent rotation of the rf field during the cavity 
passage. All fields are shown for t = 0. The cross-sections on the right show the projections in 
the x'y plane defined by a C-field tilt towards —x. The arrows illustrate the local direction (as 
given by /?) of the cavity field, but not its amplitude. As described in the text, there is a clockwise 
rotation of the field direction for an upwards passage at y < 0, and a counter-clockwise rotation at 
y>0. 

Note that the origin of the coordinate system is at the center-point of the cavity and the arrows 
shown just indicate axis orientation. 

for the subset of atoms contributing to frequency measurements. 



A. Calculating transition probability spectra 

For low transition probabilities, the Bloch equations can be solved approximately by 
assuming a non-depleted initial state population (see appendix |A]). The resulting transition 
probability spectrum is then given by 



(5) 



where f2[Aa;] is the Fourier transform of the applied pulse given in terms of the complex 
Rabi frequency Q[t] = \£l[t] \ e~ llfi , where (p describes the phase of the effective rf field. The 
coefficient q allows us to handle the variation of transition strengths [15]: For simplicity, we 
will express all field amplitudes in terms of the Rabi frequency for the [0 to 0] transition, 
such that Coo — 1- F° r the [0 to -1] transition driven by a linearly polarized rf field orthogonal 
to the C-field, the transition strength coefficient is coi = 5/16 16j . which already takes into 
account that only one of the circular polarization contributions has an effect. 

In the presence of a C-field tilt a (see Fig. [I]), the orthogonal pulse component takes the 
form 

r2j_ [t] = Vt x [t] cos a + Vt y [t] + Vt z [t] sin a . (6) 

The projected components of Q x and fl z are oriented towards the x'-direction and their 
phase will serve as a reference. Relative to this, the linear polarization of the Q y component 
is oriented at an angle of /3 = tt/2. The resulting negative phase shift ip = —71/ 2 for the 
[0 to -1] transition then yields an extra factor of % for the complex Rabi frequency fl y . 

For a single vertical passage, assumed to occur at constant velocity, the individual pulse 
components can now be expressed based on Eq. ([[]) and Eq. (EJ) as 
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where the rectangle function n ensures that there is no effect of the rf signal when the atom 
is outside the cavity. The common amplitude factor (ndbQ / (2t c cosa) is chosen such that 
for a resonant pulse with b = 1, ( = 1 and a = 0, an atom passing the cavity at p = 
experiences a n/2 pulse area for the [0 to 0] transition. The meaning of the amplitude factor 
b and the cavity specific constant ( is elucidated below. 

The Fourier-transforms of the pulse components can now be calculated as 
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describe the dependence on the detuning Au. Due to the linearity of the Fourier transform, 
the [0 to -1] transition probability given by Eq. (JSJ) now depends on the linear combination 
of the components according to Eq. as 
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For a = the expression reverts to the expected radial symmetry. 

We now explicate the role of the factors b and (: These describe the actual rf amplitude 
in relation to the value used in normal fountain operation, where it is optimized to maximize 
the contrast of the Ramsey fringes of the [0 to 0] transition. This transition is excited by 
the pulse component parallel to the C-field 



Q» [t] = Q z [t] cos a — Q x [t] sin a 



(16) 



As shown in appendix [Aj the resonant tipping angle after ge through the assumed rf 

field distribution in the cavity is then given by 
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ftii[0] = n z [0] cosa=6C|J [7P] 



(17) 



where no Q x term appears since U = for Au = 0. Averaged over the entire atomic cloud, 
this gives a transition probability after the first cavity passage of 



^oo = \ fo a (! " cos [Kf Hip] ] ) W r [p] dp 



(18) 



where the normalized distribution W r [p] describes the probability of finding any given atom 
at a position with the radial coordinate p. 

Assuming no correlation between the atomic positions during the first and second passage, 
a local maximum of the fringe contrast as a function of rf amplitude is achieved whenever 
Poo — 0.5. In analogy with [3], the cavity specific constant ( is now chosen such that the first 
of these maxima occurs at b = 1 for a uniform density distribution (where W r [p] = 2p / r^). 
For the cavity used in CSF2, this yields ( = 1.080. 
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FIG. 3. Variation of the [0 to -1] spectrum with C-field tilt angle a for a single atom traversing 
the cavity at a fixed horizontal position r, calculated using Eq. (|15p . Heavy black curve: At 
a = 0°, the spectrum is symmetric around resonance, with two lobes of equal amplitude near 
iw ro t = 2-7T x 47.2 Hz. Green shaded, dot-dashed curve: The horizontal position was chosen such 
that at a = 2.4° the field direction undergoes a rotation of constant angular velocity, resulting in 
a spectrum that is again symmetric, but with a single central lobe shifted by w ro t- Red, dashed 
curve: For an intermediate value of a = 1.2°, the lobes of the spectrum remain at similar positions 
as without tilt, but the transition probabilities are enhanced around w ro t- 

The density in a realistic atomic distribution is generally higher in the stronger field near 
the cavity axis, and a detailed analysis for CSF2 [6] yields a value of b = 0.97 for the first 
maximum of the contrast. The experiments presented throughout the paper were performed 
at a nominal 57r/2 pulse area (b nom = 5), for which we assume a five times larger amplitude, 
as described by b = 4.85. 

Note that the rf amplitude is adjusted for maximum contrast while the C-field tilt is 
already present. The necessary compensation for the reduced field component parallel to 
the C-field that can be seen in Eq. ffTB]) then causes the field amplitudes given by Eq. (J7J)- 
Eq. ([9]) to increase with the value of a even if b is constant. 

The [0 to -1] transition probability spectrum for any given microwave amplitude, C-field 
tilt, and horizontal position of cavity passage can now be found from Eq. (II 5p with an 
accuracy that is limited by the non-depletion approximation and the assumed analytical 
field distribution. Examples are shown in Fig. |3j 
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B. Estimating the angle of C-field tilt 



A point of particular interest in the spectrum of the [0 to -1] transition is the resonant 
transition probability Pgf* = Poi[0]- Since U — for Acu = 0, this depends only on the last 
term of Eq. (|T5|) as 



The non-depletion approximation is quite accurate here since the resonant transition prob- 
ability for the [0 to -1] transition is generally very low. 

As shown in appendix [Bj the C-field tilt angle a can then be calculated from the resonant 
transition probability Pgf s averaged over the entire cloud as 



even if the exact details of the density distribution in the cavity are not known. 

When the spectrum is symmetric, Pyf s is simply the transition probability at the central 
minimum and therefore easy to determine. This is not true for the asymmetric case, where 
the minimum appears shifted away from resonance (see Fig.[3]and Fig.HJ). While it is possible 
to find the proper resonant frequency from a non-linear fit of the spectrum, it is preferable 
to obtain it without introducing a new free parameter by measuring the resonant frequency 
of the [-1 to -1] transition and using the Breit-Rabi formula [17]. A series expansion yields 
the relevant terms 



where gj and gj are the nuclear and electronic g-factors, respectively, and the resonant 
frequencies Aa>[oto-i] and Aa;[_i to -i] are given relative to the [0 to 0] clock transition. Note 
that the magnetic field in the cavity can differ noticeably from the value averaged over the 
entire atomic trajectory that is measured in atomic fountains to correct for the quadratic 
Zeeman shift. 

C. Extracting position data from measurements 

Spectra calculated using Eq. ( IT5|) show that the asymmetry mostly depends on y. In this 
section we develop this into a method to extract the center-of-mass position along the y-axis 
of the atomic cloud in the cavity from a measured spectrum. 
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Summing real- valued and imaginary parts of Eq. ( Fl5|) in quadrature yields 
P 01 [Aco]= C -f(b() 2 
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The last term with its dependency on y is the most interesting and can be isolated by 
exploiting the different symmetries of the terms: Since U[—Au] = — U[Au] and G[— Aw] = 
G[Aw], the terms containing U 2 and G 2 are symmetric in Aou, while the final term depends on 
(U G) and is antisymmetric. The unwanted terms can be removed by defining the asymmetry 
of the transition probability spectrum as 

AP[Auo}= \ (Poif-Ao;] - P i[Aw]) (23) 



cos" 5 a 256 a { 7p J 



(24) 



The term in curly brackets approximates to 1 when the radial coordinate p is small, showing 
the asymmetry AP to be near-linear in the y-coordinate of the cavity passage. 

Any measured asymmetry value AP obtained in an atomic fountain is an average of many 
individual atomic spectra in the same way as in Eq. (fl~8l . except that the assumption of 
radial symmetry is no longer valid here: 

AP = J J W[x,y) AP[x,y) dxdy . (25) 

A ac 

The normalized density distribution W[x, y] is defined to be zero outside the accessible area 
A ac , given by the radius r a of the cutoff tubes. 

If AP were entirely linear in y and independent of x, this could be rewritten as the 
product of a constant slope factor s and the center-of-mass coordinate y for the ensemble of 
detected atoms: 

AP = s xy with y— / / W[x, y] y dxdy . (26) 
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The position y could then easily be determined by taking the asymmetry of a measured 
single-passage spectrum and dividing it by the value calculated for s. 

This approach does in fact provide useful results as experimentally confirmed in sec- 
tion IIVBI and discussed further in section IV Al The deviation of AP from linearity is quite 
small, and its local variation is partially averaged out since the atomic cloud generally has 
a smooth density distribution throughout the cavity, such that the variation of AP with y 
can be described through an effective slope value s e ff. 

Measuring y is then possible with the following procedure: A C-field tilt along the x-axis 
is applied. The fountain is configured such that a microwave pulse is applied only during 
a single cavity passage and ideally the microwave amplitude is reoptimized accordingly. 
The transition probability Poi is measured on resonance and at two symmetric detunings 
Aui = ±Aa; m . The resonant transition probability P^f 3 yields a value for the C-field tilt 
a Eq. (120]) . which is used to find the slope factor s e s- Finally, the center-of-mass y-position 
is extracted from the asymmetry as y = AP/s c g. The coordinate x can be obtained in the 
same way by applying a C-field tilt along the y-axis. 

IV. EXPERIMENTAL TEST 
A. Shape of observed spectra 

Even without knowledge of the exact value of s e g, the measured spectral asymmetry can 
be used to adjust the position of the launched cloud in the cavity. We used this method for 
preliminary centering of the cloud while loading atoms from a slow caesium beam. Since 
there is no centering force in the optical molasses employed in CSF2, the unidirectional 
loading from the beam tends to create an asymmetric, off-center cloud. By investigating 
the asymmetry of the spectra obtained with a microwave signal applied during either the 
upwards or the downwards passage, both the initial position of the cloud and the launch 
direction can be adjusted in this way. 

Fig. H] shows spectra taken during the centering process, one for a far off-center position 
and one after adjustment. The microwave pulse was applied during the upwards passage, 
with an amplitude that corresponds to b nom = 5. The C-field tilt was calculated as a = 2.4° 
from the resonant transition probability. Note that the C-field tilt inside the cavity of CSF2 
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FIG. 4. Comparison of experimental (symbols) and calculated single-passage spectra (lines). 
Experimental data was taken with the initial cloud deliberately placed off-center (red circles) and 
with the molasses beam balance adjusted for improved centering (green triangles). This data was 
then fitted with Eq. (|15p by varying the horizontal position of passage r for a C-field tilt of 2.4° 
extracted from the resonant transition probability (red, green lines). To show the effect of the 
C-field tilt, the spectrum for centered adjustment was recalculated assuming a = (black line). 
Two choices of the measurement detuning (see section IIV Ap have been marked. Dashed vertical 
lines indicate oo m = ±w rot (maximum sensitivity) and solid lines indicate u m = ±4.3/£ c (suppressed 
effect of frequency uncertainty). 

so far cannot be deliberately controlled and simply results from the superposition of stray 
magnetic fields with the homogeneous, vertical field produced by the C-field coil. 

To demonstrate the adequacy of the analytical model, the figure also shows a non-linear 
fit of each spectrum with Eq. (1151) . which describes the spectrum for a single atom. The 
best agreement is achieved for a cavity passage at a horizontal position with y oS = 1.9 mm 
in the off-center case and y ccnt = 0.0 mm after centering. 

The highest sensitivity for a position measurement is obtained by measuring at a detuning 
Au m = Wrot, where U G = —1, maximizing the value of AP. However, as shown by the 
dashed vertical lines in Fig. HI the transition probability is then measured in a region of the 
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spectrum with a considerable slope. This translates any error in the determined resonant 
frequency into an error of AP that can be large: In the off-center example shown in the figure, 
a realistic frequency uncertainty of <7 res = 2tt x 1 Hz causes an uncertainty contribution of 
capm = 1-5 x 10~ 3 , which is considerably larger than the statistical measurement uncertainty 
A good choice for an optimized measurement detuning is Au m = 4.3/t c , such that the 
measurements are effectively taken at the maxima of U (see Eq. (fT5|) ). which provides the 
largest contributions to the transition probability spectrum. Using this for the sample 
measurement reduces the uncertainty contribution to an insignificant value, at the price of 
reducing the position sensitivity by 13%. 

B. Variation of AP with launch direction 

To further confirm the calculations, a series of experiments was performed. One of the 
design features of the fountain CSF2 is a goniometric stage for the entire molasses zone 
that makes it possible to vary the launch direction in a highly reproducible fashion without 
tilting the entire fountain 

Measurements of the detected atom number, the resonant transition probability PJf 3 
and the asymmetry AP were taken for a range of launch angles up to 3 mrad, covering 
both the position sensitive ^-direction and the insensitive x-direction. Since it is unknown 
which alignment results in a vertical launch, the tilt of the launch direction is given by the 
angles S x and S y relative to the direction that maximizes the returned atom number. This 
optimized tilt could be determined to within an uncertainty of as = 0.05 mrad by fitting 
the atom number detected as a function of the launch tilt in x- and y-direction with a 2D 
Gaussian distribution for a suitable subset of measurements. The measurement detuning 
for the determination of AP was chosen as Aoo m = 2ir x 65.2 Hz, based on the optimized 
value from section HV Al All measurements were performed twice, with the microwave pulse 
applied either during the upwards passage or during the downwards passage, and the results 
given are averaged over 100 successive launches. 

After correcting for crosstalk in the detection system, the measured resonant transition 
probabilities Pq^ s for the upwards and the downwards passage at each launch tilt all agree 
with the average value of P™ s g = 2.86 x 10~ 3 to within a measurement uncertainty of op = 
5 x 10~ 4 . For a nominal pulse area described by b nom = 5 and using Eq. (120"]). this result 
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corresponds to a tilt angle of a = 1.4°. The difference to the tilt angle of 2.4° found in 
experiments conducted several months earlier (such as the one in section IIV Aj) indicates a 



change in the stray magnetic fields that is likely due to work on the attached caesium beam 
apparatus over the intervening period. 

Looking at the spectral asymmetries AP[5 X , 5 y ] for either the upwards or the downwards 
cavity passage, a three-dimensional plot of the data points obtained shows them to lie within 
a plane. As the tilt of the C-field used in the measurements was originally unintended and of 
unknown direction, the orientation of the plane fitted to the data for the downwards passage 
was used to align the coordinate system such that a change of the launch tilt S y (along the 
2/-axis) results in the largest change of AP. 

Fig. G2a) can be regarded as a cross-section of the three-dimensional plot: It shows AP 
for all investigated launch angles plotted over the relative launch tilt 5 y . The deviations from 
the linear fit are then independent of 5 X and consistent with a purely statistical variation 
with standard deviations of cr^p = 4.3 x 10~ 4 and a^p = 3.6 x 10~ 4 for the upwards and 
downwards passages. The slight increase over the value expected from simple uncertainty 
propagation, o"ap = ap/y/2 « 3.5 x 10~ 4 , is attributed to the decrease in detected atom 
number at large launch tilts and the resulting reduction in signal-to-noise ratio. 

Using the effective slope values s^ ff = — 6.35m~ 1 and s^ ff = 6.26 m _1 (obtained in sec- 
tion |IVC]) now yields a linear fit for the center-of-mass positions y as a function of <L: 



y*[6y] = -0.47 mm - (0.14 mm/mrad) S y (27) 
y±[S y ] = -0.01 mm + (0.28 mm/mrad) 5 y (28) 

for the upwards and downwards cavity passage, respectively. 

The fact that y^[5 y ] < for all tested launch tilts indicates that the cloud in the molasses 
zone is positioned off-center even before launch. Due to the lack of a centering force, this 
can easily happen if the cooling beams are not perfectly aligned. Adjusting the launch 
tilt then mostly affects the cloud position during the downwards passage 2r[5y], which the 
measurements confirm to be well-centered at 5 y = 0, where the maximum returned atom 
number is detected. 

Interestingly, the cloud position on the way up shifts against the launch tilt due to a 
position dependent selection process: For a launch tilt 5 y > 0, atoms with an initial position 
y > are more likely to drift out of the accessible area before the second cavity passage and 
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FIG. 5. (a) Measured asymmetry values AP over launch tilt 5 y in y-direction together with linear 
fits. Green upward triangles mark data taken for the upwards passage, blue downward triangles for 
the downwards passage. The red dashed line marks the tilt setting for a vertical launch direction 
obtained from the simulation, (b) Combinations of launch tilts 5 X and 5 y where data was taken. 
The red x marks the direction of vertical launch. 



are therefore less likely to be detected. Conversely, atoms placed at y < now tend to move 
towards the fountain axis which increases their chance of detection. For the ensemble of 
atoms that eventually get detected, the center-of-mass position during the upwards cavity 
passage then appears shifted towards — y, against the launch tilt (see also section 6 of Q]). 
This shift is not compensated by the additional horizontal motion between launch and first 
cavity passage unless the cloud is very small. 

The same selection process also causes a general suppression of the cloud shift with tilt: 
For 5 y > 0, atoms with an opposite thermal velocity component t>* h < are more likely 
to reach the detection zone than those with -y* h > and therefore contribute more to the 
detected ensemble. This is the reason why the observed shift of the cloud between passages 
of 0.42 mm per mrad launch tilt is much smaller than the value of 2.5 mm per mrad expected 
for a simple parabolic trajectory. 
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C. Numerical simulation 



To confirm the effect of these selection mechanisms on the observed center-of-mass posi- 
tions, a numerical trajectory simulation was developed. It is based on Monte-Carlo methods 
and generates atomic trajectories starting from a three-dimensional cloud with a Gaussian 
distribution of initial positions and velocities. Trajectories are rejected if they intersect any 
of the known restrictions, most importantly the cutoff tubes of state selection and Ramsey 
cavities, or if they miss the detection zone. 

For each remaining trajectory, the horizontal position of the upwards or downwards cav- 
ity passages is determined, such that AP for the upwards and downwards passage can be 
obtained from the average transition probabilities at the detunings ±Auj m . To eliminate the 
error introduced by the non-depletion approximation, the transition probabilities are found 
by numerically solving the differential equations Eq. ( ]A2j) and Eq. (1A3j) instead of using 
Eq. (H5J). 

The best agreement with the experimental data for AP and the detected atom number, 
both as a function of launch tilt, is achieved for an initial cloud with a horizontal density dis- 
tribution described by a T = 3.0 mm, placed at a horizontal position offset of yo = —1.1 mm. 
The horizontal velocity distribution is given by a v = 7.9mm/s, corresponding to a temper- 
ature of T c i ou d = 1.0 /iK. 

Even at this very low cloud spread and for optimized launch direction, less than 15% 
of the launched atoms are detected at the end of the cycle. Any change in the launch tilt 
then strongly affects which fraction of the cloud is most likely to be detected, confirming 
the counter-intuitive shift of the cloud on the upwards cavity passage and the reduced shift 
on the downwards passage. 

After optimization, the trajectory simulation also provides values for the effective slope 
value s e fj (compare section IIII Ch : For a range of launch tilts varying in both x- and y- 
direction, we calculate the center-of-mass position y of the simulated atomic ensemble in 
addition to AP. By plotting AP over y, the effective values s^ s = —6.35 m -1 and s^ s = 
6.26 m _1 for the upwards and downwards passage are then obtained from linear fits. The 
sign reversal for s^ s results from the vertical antisymmetry of H p . 
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D. Assumed density distribution 

The Monte-Carlo nature of the trajectory simulation complicates a systematic evaluation 
of the position uncertainty arising from the distribution of atoms in the cavity. Because of 
this, we introduce a simplified model to describe the density distribution inside the cavity 
during the upwards or downwards passage, for those atoms that are detected at the end of 
the cycle and are therefore relevant for frequency measurements in the fountain. 

The trajectory simulation shows that due to clipping at the cutoff tubes and the limits 
of the detection zone, the density of these relevant atoms falls off sharply in the region near 
the edge of the accessible area, so that we can describe the relevant atomic density by the 
approximate distribution 



which is constrained to a radius r r, m around the cavity axis. The density fall-off within this 
radius (which is smaller than the cutoff tube radius r a ) is described by a, while x c and y c 
provide a way to vary the center-of-mass position of the cloud. The normalization factor a n 
is calculated numerically. 

Values for AP and y and can now be found from Eq. (12 5 j) and Eq. (1261) . In the former, 
it is again preferable to calculate the local values of AP from a numerical solution of the 
differential equations. Since both y and AP are zero for a centered cloud (where x c = y c = 0), 
the effective slope value for an exemplary deflection of the cloud is 



The exact values chosen for x c and y c are not critical, as the value obtained for s opt changes 
by less than 0.5% for center-of-mass positions within 1mm of the origin. In the following, 
we will assume s opt to be calculated for a near-zero deflection. 

The approximate distribution W a reproduces the slope values obtained directly from the 
trajectory simulation for a common limiting radius ry, m = 4.0 mm, if "undipped" cloud sizes 
of = 3.3 mm and cr^ = 7.3 mm are chosen for the upwards and downwards passage. These 
cloud sizes are in good agreement with the expected thermal expansion of the launched 
cloud. 




(29) 



■^opt 



AP[s c ,y c ] 

y[x c , Vc] 



(30) 
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V. UNCERTAINTY ESTIMATE 



Investigating the dependence of the spectral asymmetry on the cloud position was made 
possible by an uncontrolled and originally unintended tilt of the C-field in CSF2. This does 
not provide position information along the currently insensitive x-axis, and the accuracy of 
the position measurements along the y-axis is also not optimal. 

For the discussion of the achievable uncertainties, we will therefore assume a hypothetical 
version of CSF2 in which a set of coils provides an intentional tilt of the C-field that can 
be changed in both amount and direction. This allows operation at optimized parameters, 
which will be set as follows throughout the section: Choosing an even value for b nom makes it 
possible to adjust the rf amplitude for maximum [0 to 0] transition probability after a single 
cavity passage, as the tipping angle is then a multiple of tt. Compared to the optimization 
after both cavity passages that otherwise needs to be used, this agrees more closely with the 
assumptions used in the derivation of Eq. (1201) and improves the accuracy of the determined 
C-field tilt a. It is preferable to reach an even 6 nom through a reduction of the nominal pulse 
area to 6 nom = 4 (with b ps 3.9), since this also reduces the overall transition probability, 
which improves the linearity of AP in y. An increase of the C-field tilt to a value close 
to a = 3° more than compensates the reduction of the slope value due to the reduced rf 
amplitude to give s op t ~ 8.9 m"" 1 . 

A. Residual position uncertainty after centering 

Ideally, the asymmetry measurements are used to center the atomic cloud in the cavity 
during both the upwards and downwards passage, and in both x- and y-directions. When 
AP is zero, the uncertainty of the calculated slope factor becomes insignificant. If the 
observed spectra resulted from a single atom, the residual uncertainty would then be a S y ng = 
| s~pt °ap| = 0.04 mm. 

Reducing cap through increased averaging is only possible to the point where systematic 
effects such as the nonlinearity of the detection system or the transition probability back- 
ground from the off-resonant excitation of the [0 to 0] transition become limiting. For the 
parameters used here, the latter was estimated to cause an error in the measured AP value 
on the order of 2 x 10 -4 , based on numerical calculations. At a similar level, the uncertainty 
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FIG. 6. Variation of the spectral asymmetry with the horizontal position of cavity passage given 
by x and y. Values of AP[x,y] are given by the contour labels for optimized parameters a = 3.0°, 
b = 3.9 and Aw m = 2tt x 65.2 Hz. The variation of AP with x can cause an error in the position 
determination from AP that is at its maximum for the situation illustrated by the marked points 
(see text). 

resulting from a possible error in the assumed resonance frequency of the [0 to -1] transition 
discussed in section [TV Al might also need to be included in <7y ing . 

For an extended distribution of atoms there is another, larger contribution to the un- 
certainty: A value of AP = is the result of positive and negative AP contributions from 
individual atoms. Since AP is not entirely linear in y and also has a dependence on x, 
a density distribution with a center-of-mass position of y = does not necessarily yield 
AP = and vice versa. 

Fig. [6] illustrates this: The contour plot shows the variation of AP with the horizontal 
position of cavity passage for a single atom. Clearly, the same value of AP can occur 
for a range of y-positions if x is varied. To quantify the resulting position uncertainty, 
we investigate a minimalistic distribution of just two atoms placed such that y = and 
therefore y\ = —y%. The largest deviation of AP from zero occurs when the atoms are 
placed as marked in the figure, one at x± — and the other at the edge of the cutoff tube. 
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The resulting asymmetry of AP = —2.7 x corresponds to an error of the extracted 
center-of-mass position of Ay = 0.31mm. 

This two-position distribution is the worst-case arrangement and allowing other distribu- 
tions with more atoms will not increase the resulting error. Fortunately the atomic density 
distribution normally is not concentrated in localized regions, but spread smoothly across the 
entire accessible area. Simply requiring a constant density distribution along the connect- 
ing line between the worst-case positions, as indicated in the figure, reduces the maximum 
position error to Ay = 0.11mm, which we will use as the distribution based contribution 
to the uncertainty of the center-of-mass position after centering, (jy lst . Together with cr^ ing , 
this yields a combined centering uncertainty of 0.12 mm for a measurement at optimized 
parameters. 



B. Position uncertainty for an uncentered cloud 

We now investigate how critical the assumed values for cr , and rn m are. Since the 
center-of-mass position is extracted from the measured asymmetry as y = AP/s, we look at 
the uncertainty of the inverse slope factor s~p t . 

Simultaneously varying cfi (or cH") by 25% and the limiting radius r\ im from 3 mm to the 
full cutoff-tube radius of 5 mm changes the calculated value of s~\ by a maximum of 8.5%. 
We will use this number as a conservative estimate of the uncertainty of the inverse slope 
factor due to the density distribution of atoms in the cavity. This result once again assumes 
a numerical solution of the differential equations. Calculating the effective slope values from 
Eq. (1231 directly would introduce an additional error of around 6% due to the non-depletion 
approximation. 

A comparable uncertainty contribution arises due to the use of the approximated analytic 
field distribution of Eq. ([Q) and Eq. (T5]): The radial variation of the tipping angle for the 
actual fields in a cavity with endcap apertures is (sj 

e act [p] = ^&Jo[M , (3i) 

where k = u/c is the wave vector of the applied frequency. The factor rj is a constant 
that depends only on the cutoff tube radius and is equivalent to ( in Eq. ©-Eq. ([9]). 
For CSF2, the expression for O derived from the approximate fields and given in Eq. ( TT7|) 
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underestimates the tipping angle by 3.6% on the cavity axis and overestimates it by 4.8% 
at the cutoff tube radius. This provides a measure of the components of the actual H z fields 
that arise due to the presence of the endcap apertures and are neglected in the analytical 
model. For H p , the neglected components can be assumed to be of similar magnitude. Since 
the asymmetry measurements are taken near the extrema of the Fourier transform H p (see 
section IIV A|) . it is unlikely that the neglected field components have over-proportionally 
large Fourier components at these frequencies. 

Realizing that AP is essentially given by the product of H p and H z (which provide the 
U and G terms of Eq. f[T5"j) ) then gives a maximum expected increase of the asymmetry by 
(1.048) 2 . We will represent this as an assumed 10% uncertainty of the inverse slope factor. 
The value depends on the aspect ratio of the cavity, however: While the 7 term in Eq. ( TTT|) 
varies with the cavity radius, k in Eq. f l3Tj) does not. 

The final significant uncertainty contribution results from the uncertainty of the param- 
eters used to calculate the effective slope factor. This can be estimated from Eq. ( 1241) : 
Approximating the term in brackets as 1 and combining the result with Eq. (1201) . this yields 



Standard uncertainty propagation can now be applied. For simplicity we assume cos a = 1 
and parameter values typical for CSF2. A predicted value for Pq° s at optimized parameters is 
calculated from Eq. ffT9]) . Significant uncertainty contributions arise from the measurement 
of P, with an assumed uncertainty of op = 5 x 10 -4 (see section IIVBI) . and from the 
uncertainty of the rf amplitude <7b = 0.1. In total, the calculation yields a relative uncertainty 
of the inverse slope factor of 6.4%. 

Adding all three contributions in quadrature, we assume a total relative uncertainty of 
u s = 15% for the inverse slope factor at optimized parameters. 

Similar arguments as in section IV Al can be made for the part of the position uncertainty 
that does not depend on the slope factor. Assuming the same values for cr® ing and a^ mt as 
before, the total uncertainty is then given by 




(32) 




(33) 



This yields cr* ot = 0.14 mm for a cloud with a center-of-mass position of y = 0.5 mm, as 
encountered during the upwards passage in CSF2. 
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VI. SUMMARY 



We have demonstrated a method to obtain information on the position of atoms in the 
cavity of a fountain clock by applying a slight tilt of the normally vertical C-field. Since 
measurements are performed using the actual field in the cavity and the actual detection 
system, the results represent the position of the subset of atoms that contribute to frequency 
measurements. 

Although the C-field tilt in the fountain CSF2 was originally unintended, the applica- 
tion of this method provided valuable information on the cloud positions during the two 
cavity passages, which has already been used in a new DCP evaluation that was recently 
published jfj]. 

Systematic position measurements will be possible for any fountain clock that is equipped 
with a means to provide a controlled horizontal component of the magnetic field on the order 
of 10 nT in the cavity region. The methods shown here then allow for efficient measurements 
of the cloud position that might even be periodically performed between frequency measure- 
ments. Since the results also provide information on the precise value of the C-field tilt, they 
can also be used to exclude an unwanted tilt during the actual frequency measurements. 
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Appendix A: Quantum mechanical derivation of the transition probability 

Following a standard textbook approach, we start from a two-level model with the eigen- 
states ) an d hAi)? such that the state as a function of time can be written as the linear 
combination 



The coefficients Co and c\ are generally complex and incorporate information on the phase 
of the atomic state. The behavior of the system in the presence of an oscillating field can, 
after applying the rotating wave approximation, be described by the two coupled differential 




(Al) 
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equations 

c' [t] = -i 1 0* [t] exp [i Au t] c x [t] ( A2) 

c[[t] = -i\ Q[t] exp [-i Aw t] co[t] . (A3) 

Here Aw = u — u ms is the detuning of the microwave signal from the atomic resonance. 
The rf field amplitude and phase (p are described in terms of the complex Rabi frequency Q 
(see section Ull Ap . An analytic solution exists for the case of a resonant pulse with constant 
phase: For an initial state described by Co[tj] = 1 and c\\tj\ = 0, the transition probability 
depends on 



Q[t] dt 



(A4) 



Ci[tf] = —i sin|6 , with = 

where will be referred to here as the resonant tipping angle. The probability P of finding 
the atom in the excited state at tf, after the end of the pulse is then 

p= \ Cl [t f ]\ 2 = (1 -cos©)/2 . (A5) 

Insight into the off-resonant transition probabilities can be gained by using a non- 
depletion approximation: For small transition probabilities Cq is approximately 1, so that 
Eq. (1A3j) can be integrated to give 



cf a \t f ] = -i\ \ Q [t] exp [-% Aoo t] dt . ( A6) 

By letting ti and tf go to — oo and oo respectively, the integral becomes identical to the 
Fourier transform of Q[t], where we use the convention 

/oo 
Q[t] exp[-z Aut] dt , (A7) 
-oo 

because it directly includes the resonant tipping angle as = Q[0] , which is valid even 
outside the non-depletion approximation. 

Appendix B: Derivation of the tilt angle equation 

Averaged over the radial distribution W r [p] introduced in Eq. ( ITS]) , the resonant [0 to -1] 
transition probability given by Eq. f[T^]) becomes 

PoT = W ^ (&Cf MlP] tan«) 2 H/ r [p] dp . (Bl) 
Jo 
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Using (Jq[x\ ) 2 = 2J [x] — 1 + 0[x] 4 , this yields 

% s = \ coi (K f tan a) 2 ( 2 ° J [7P] Wr [p] dp 



£w r [p]dp\, (B2) 



where the small fourth order term has been ignored. The normalization of W r [p] makes the 
second integral 1. And since the same distribution was used to optimize the rf amplitude, 
the approximation 

f / Q\p]W r \p] dp 

Jo 

bClUlp] W r [p] dp (B3) 
can be used to replace the first integral and find 

= J cm (f tana) 2 (6 nom + A6) (b nom - A6) . (B4) 
Since A6 = b( — b nom is small, this gives 

% S = I c oi(f fenomtana) 2 . (B5) 
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